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Abstract. The large complex zeros of the Jost function (poles of the S matrix) 
in the complex wave number-plane for s-wave scattering by truncated potentials are 
associated to the distribution of large prime numbers {p n } as well as to the asymptotic 
behavior of the imaginary parts {t n } of the zeros of the Riemann zeta function on 
the critical line. A variant of the Hilbert and Polya conjecture is proposed and 
considerations about the Dirac sea as "virtual resonances" are briefly discussed. 
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^ '. PACS numbers: 02.10.Lh, 03.65.Nk, 11.55.Bq 

There is an old conjecture attributed to Hilbert and Polya about the zeros of the 
Riemann zeta function on the critical line as eigenvalues of a self-adjoint linear 
operator H in some Hilbert space. Ever since Montgomery's pQ discovery of these 
zeros behaving like the eigenvalues of a random hermitian matrix, attempts have been 
made in order to find a quantum system with the Hamiltonian represented by such 
an operator (see Berry and Keating j2j, and references therein). We report in this 
Letter a variant for the above conjecture: instead of looking for H, whose spectrum 
coincides with the Riemann zeta zeros, we are looking for complex momenta poles of the 
scattering S matrix such that by a given transformation they are all mapped into the axis 
Re z = 1/2 in coincidence with the Riemann zeroes. The associated quantum system 
could be the "vacuum" , interpreted as an infinity of "virtual resonances" , described by 
the corresponding S matrix poles. 

The first to associate the Riemann hypothesis in terms of transient states were 
Povlov and Fadeev [3] by relating the nontrivial zeros of the zeta function to the complex 
poles of the scattering matrix of a particle on a surface of negative curvature. Here the 
same problem is discussed in the usual scattering by a potential. We begin by showing 
that complex momenta zeros of the Jost function for s-wave non-relativistic scattering by 
repulsive cutoff potentials, after appropriate transformations, correspond to the global 
behavior of heights {t n } of the zeros of ((z) on the critical line, z n = 1/2 + i t n , and 
surprisingly to the global behavior of the prime sequence {p n }- With the aid of these 
transformations it is possible to obtain an approximate formula connecting the nth 
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prime and the nth zeta zero, which we believe will be useful to strategies for primality. 
The local behavior, defined as deviations from the average density of the zeta zeros, are 
not obtained by the potential considered here. This will be the object of a forthcoming 
work where a statistical hypothesis with respect to some residual interaction will be 
introduced jlj. Actually fluctuations in {t n } are interesting for their universality, being 
observed in quantal spectra in different physical systems (see Mehta jS]), and by the 
connection with chaotic dynamics, (for a review see Bohigas [6 ). 

The distribution of primes, {p n }, among all natural numbers, n, in spite of their 
local deviation of any known order when viewed at large possesses regularities that can 
be approximated by some formulas. The approximate number of primes tt(x) less than a 
given x, also called the prime counting function, is given by the prime number theorem 
7r(x) ~ xj lnx (see Titchmarsh Chapter- III), where lnx is the natural logarithm of 
x. This relation gives the asymptotic approximation for nth prime p n , 

p n ~ n In n , as n — > oo . ( 1 ) 

The connection between the distribution of prime numbers ir(x) and complex zeros 
of the zeta function, a kind of duality between the continuum and discrete in number 
theory, started with Riemann's 1859 paper (see Edwards jH], p. 299) by introducing 
methods of analytic function into number theory. Riemann's zeta function is defined 
(Ref. |7j, p.l) either by the Dirichlet series or by the Euler product 

((z) = J2n- z = U( 1 -P~T 1 , Rez>l, (2) 

n p 

where n runs through all integers and p runs over all primes, ((z) can be analytically 
continued to the whole complex plane, except at z — 1 where it has a simple pole with 
residue 1. It satisfies the functional equation ((z) = 2 z tt z ~ 1 sin(7T2:/2)r(l — z)C(l — z) , 
called the reflection formula, where V{z) is the gamma function. It is known that ((z) 
has simple zeros at points z = —2n, n = 1, 2, . . ., which are called trivial zeros, with an 
infinity of complex zeros lying in the strip < Re z < 1. From the reflection formula 
they are symmetrically situated with respect to axis Re z = 1/2, and since ((z*) = C*( z ), 
they are also symmetric about the real axis, so, it suffices to consider the zeros in the 
upper half of strip 1/2 < Re z < 1. It is possible to enumerate these complex zeros as 
z n = s n + i t n , with t\ < t 2 < t 3 < . . ., and the following result can be proven (see 
Titchmarsh 7^ , p. 214) 

i 27m 

\z n \ ~ t n ~ , as n > oo. (3) 

Inn 

The Riemann hypothesis is the conjecture, not yet proven, that all complex zeros of 
C(z) lie on axis Re z = 1/2, called the "critical line". Based on this conjecture Riemann 
improved on Gauss's suggestion that n(x) approximate the logarithmic integral as 
x — > oo with a new prime number formula, taking into account local prime fluctuations 
in terms of nontrivial zeta zeros (Ref.jH], p. 299). 

The Jost function |Hj has played a central role in the development of the analytic 
properties of the scattering amplitudes. In order to recall its properties, let us consider 
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the scattering of a non-relativistic particle, without spin, of mass m by a spherically 
symmetric local potential, V(r), everywhere finite, behaving at infinity as 

V(r) = 0(r- 1_e ), e > 0, r -> oo . (4) 

The Jost functions f±(k) are defined (see Newton p. 341) as the Wronskian 

W, f±(k) = W[f±(k,r), ip(k,r)], where <p(k,r) is the regular solution of the radial 
Schrodinger equation 

[^- 2 +k 2 -V(r)] ¥ >(k,r)=0, (5) 

(in units for which h = 2m — 1) k being the wave number and the Jost solutions, 
f±(k,r), are two linearly independent solutions of equation (jSJ). They satisfy the 
boundary conditions lim r _, 00 [e =R/cr /±(A;, r)] = 1, corresponding to incoming and outgoing 
waves of unit amplitude. 

The properties of the solution of differential equation © define the domain of 
analyticity of the Jost functions f±(k) on the complex fc-plane as well as its symmetry 
properties, such as for real potentials, f*(k*) = f-(k). The phase of the Jost function 
is just minus the scattering phase shift 8(k), that is f±(k) = |/±(A;)| e T%5 ^ , so that 
the usual S matrix is given by 

S(k) es e«W = . ( 6 ) 

The complex poles (Re k ^ 0) of S(k), or zeros of f+(k), correspond to the solutions of 
the Schrodinger equation with purely outgoing, or incoming, wave boundary conditions. 
Resonances show up as complex poles with negative imaginary parts, their complex 
energies being 

kn = — i (7) 

where S„ and T n represent the energy and the width, respectively, associated with nth 
resonance state. For small T n , resonances appear as long-lived quasistationary states 
populated in the scattering process. If the width is sufficiently broad no resonance effect 
will be observed, as if the lifetime of state 1/T n is smaller than the time spent by the 
particle to traverse the potential: we will call this kind of S matrix pole as "virtual 
resonances" throughout. "Virtual resonances" like the long lived observed ones, in fact, 
are represented by pairs of symmetrical S matrix poles in the complex k-plane, a capture 
state pole in the third quadrant and the decaying state in the fourth quadrant, since they 
give to the asymptotic solution an incoming growing wave and an outgoing decaying 
wave, respectively, exponential in time [TT]. That could be described by Gamow vectors 
treated by Bohm and Gadella [12] as pairs of S matrix poles corresponding to decay and 
growth states. Examples of broad resonances are the well known large poles, related 
basically to the cutoff in the potential considered to be without any physical significance 
(see Nussenzveig (THj, p. 178). 

Condition which together with differential equation (jSJ) establishes domains of 
analyticity for f+(k), is not sufficient to determine the asymptotic behavior for its zeros; 
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for that we would need more information about the interaction. Then the potential is 
set equal to zero for r > R > 0, which is the cutoff of the potential at arbitrarily 
large distances R. With this restriction it can be shown that f+(k) = is an entire 
equation of order 1/2 , and according to Piccard's theorem, has infinitely many roots for 
arbitrary values of the potential. Position k n of these zeros determines f+(k) uniquely 
in the whole complex plane, as a consequence of Hadmard's theorem, which provides 
[T4] f±(k) = e ±tkR f(0) rr?£Li(l — k/k n ). The asymptotic expansion of k n , for large n, 
after introducing dimensionless parameter {3 = kR, is given by (Ref.|l()j. p. 362) 

. (a + 2) In \n\ ,„ s 
Pn = nn-i± ^^^ + 0(1) (8) 

where n = ±1, ±2, ±3 • • • and a is to be defined by the first term of the potential 
asymptotic expansion, near r = R, through V(r) = C(R — r) c H — a > and r < R. 

The connection between the complex zeros of the Jost function and those of the 
Riemann zeta function is obtained by means of the transformation: 

by which the lower half of complex /3-plane (Re f3 ^ 0) is mapped onto the critical 
axis, Re z = 1/2, of complex z-plane. This suggests a variant for the Hilbert and 
Polya conjecture, namely, looking for a potential that gives a Jost function with all 
zeros on the lower half of complex /3-plane (Re (5 ^ 0) that coincide with complex zeros 
of the Riemannn zeta function after the transformation (JHJ). In this way the Riemann 
hypothesis follows. For real potentials, these complex f3 zeros are located symmetrically 
about the imaginary axis, then by Q they will be mapped symmetrically about the 
real axis into the critical line. 

Now we show that for cutoff potentials transformation Q gives rise to complex Jost 
zeros with the same asymptotic behavior as the complex Riemann zeta zeros, being all 
in the critical line. If {f3 n } are the zeros of f+(0) then by we get 
1 . Re (3 n 

Zn= 2- l 2J^K^ (10) 
from (jHJ), with a = 0, we see that {Im z n } has the same asymptotic expansion as {£„}, 
given by (jHJ), i.e., 

Im z„ = ^ as rnoo, (11) 

which means that for each resonance, ratio {2Re (3 n /Im (5 n } corresponds to the height 
of the zeta zero on the critical line. 

On the other hand, after introducing dimensionless quantities, energy E n = R 2 S n 
and widths G n = R 2 T n , equation ((7j) is written as 0^ = E n — iG n /2, the dimentionless 
widths {G n }, defined as 

G n = 4 Rep n Imp n , (12) 
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after taking into account (|SJ), when a = 0, shows the same asymptotic expansion for 
large primes (0), given by the prime number theorem, 

G n = Airp n as n — ► oo. (13) 

Then nth large complex Jost zeros are also related to nth large primes, showing an 
asymptotic connection between primes and complex Riemann zeta zeros, in a one-to- 
one correspondence. 

In the scattering on a surface of constant negative curvature |3| the potential is 
replaced by imposing that the particle move on the given surface in order to obtain 
the scattering function S(k) in terms of the Riemann zeta function. Provided the 
Riemann hypothesis is true the poles of S-function in the complex /c-plane are given 
by k n = t n /2 — i/A |3J EH]. By using transformation we obtain z n = 1/2 + i t n 
in coincidence with the complex zeros of the zeta function as expected. It follows also 
from the usual resonance width definition (j2J) that T n = t n /2, is therefore not constant 
as considered in jTHj, exhibiting the same fluctuation of the height of the zeta zero in 
the critical line, in accordance with the Wardlaw and Jaworski ^Hj results on the time 
delay fluctuation for this kind of unusual scattering. 

In order to look for physical systems associated to the proposed S'-function, we will 
return to the conventional scattering. In this connection it is interesting to examine 
more explicitly the potential parameter dependence in the above asymptotic expression 
(jHJ. Khuri .17] has recently proposed a modification to the inverse scattering problem 
in order to obtain the potential whose coupling constant spectrum coincides with the 
Riemann zeta zeros. The model we have chosen is the non-relativistic s-wave scattering 
by a spherically symmetric barrier potential, V(r) = Vq for r < R and for otherwise. 
From the stationary scattering solution with this potential one obtains the Jost function 
f+(k) = e tkR [k'cos(k'R) — ik sin(k'R)] , where k' 2 — Vq — k 2 . Introducing dimensionless 
parameters: a = k'R, j3 = kR and v = VqR 2 , the zeros of the Jost function are given 
by the solution of the complex transcendental equation 

\jfi 2 -v cot yflP-v = i/3 (14) 

for each value of potential strength v. The displacements of the roots, (3 m in the (3- 
plane with the variation of the potential strength are shown by Nussenzveig [15J. These 
zeros are all in the lower half of complex /?-plane and located symmetrically about the 
imaginary axis. For finite v ^ 0, in equation (|14J) . the large values of the (3 n zeros are 
asymptotically determined [HI [T^j , as 

ln(2 n-K/\jv\) 1 

f3 n = ±nn =F — - i In 2 utt U/v ) t — ) , (15) 

nix v mi 

n being a large positive integer. This formula gives the dependence of the potential 

strength v, considered as a free parameter, in the determination of the asymptotic 

Jost zeros (3 n to be compared, after the transformation (fTUjl . with the zeta heights 

t n on the critical line, and after ()12j) with asymptotic prime numbers p n . By using 

asymptotic expansion (|T5|) we calculate the energy/width ratios {AnE n /G n }, where 
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Figure 1. (a) The height of nontrivial zeta zeros on the critical line. The full line 
according to energy/width ratios {4 TiE n /G n } from the Jost zeros; the dashed line 
from the asymptotic formula ©; the points are computed by Odlyzko [18]. (b) Primes 
in order of size. The full line according to dimensionless width {G n /47r} from the Jost 
zeros; the dashed line according to the prime number theorem £Q); the dotted line 
taken from the table of Caldwell [19]. Insets are defined in the text. 

4irE n /G n ~ 7tRe [3 n /Im /3„, which are compared to {t n } computed by Odlyzko jTS] 
and shown in figure 1(a), for v = 2. Then we see the approximate agreement from the 
beginning, for n up to 6 x 10 4 . The same Jost zeros j3 n , for v — 2, after transformation 
(fT2j) . give dimensionless widths {G n /47r} that are compared with prime sequence {p n } 
taken from the table of Caldwell [19,, as shown in figure 1(b). The numerical factor 
1/47T was imposed to G n by the asymptotic limit (fT^j) from the prime number theorem. 
Here we see the agreement with the global behavior of the sequence of primes, from 
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the beginning in the same range. In the upper left hand corner in both figures the 
local zeta zero and prime fluctuations shown by dots are not present when obtained 
by the corresponding Jost zeros (open circles), here normalized on the 9880-th zero; the 
normalization factor corresponding to the height on the critical line is Nz = 1.071 and 
to the prime being Np = 0.978. The inset of each figure in the lower right hand corner 
shows deviation in the height of zeta zeros figure 1(a) and prime figure 1(b), obtained 
by the corresponding Jost zeros, from true values taken from tables [T%1 IT9"] in the range 
considered. 

Finally, a conjecture is made in order to associate the Dirac "virtual 
resonances". The hole theory is the interpretation of the negative energy solution of 
the relativistic single-particle Dirac equation in which the vacuum consists of all these 
negative energy states being filled with electrons, such that it could be considered from 
the beginning as a many-particle system being described by the formalism of the single- 
particle theory. According to Dirac (Ref.|22|, P-34), "The vacuum must be a state with 
a lot of particles present corresponding to some stationary solutions of the Schrodinger 
equation. But there are no known solutions of this Schrodinger equation - not even a 
solution which could represent the vacuum" . The question about the vacuum structure 
was bypassed by the second quantization, where a vacuum state is assumed and an 
operator defined, in order to create an electron when applied to this state, without 
knowing what the vacuum state really is. We would like to suggest a description of 
the vacuum structure as being a dynamical system described by "virtual resonances", 
completely independent of the second quantization. Instead of looking for a stationary 
solution, we look for a transient state of very short time duration. Specifically we propose 
a dynamical model for the vacuum described by infinite denumbered virtual resonances, 
with discrete widths and energies, which could be useful in the description of quantum 
chaos, which is under investigation. It could perhaps be verified experimentally by some 
devise that amplifies the vacuum fluctuation; the fluctuating nature of the Casimir 
force was recently discussed by Bartolo et al. [201 • I n this case, the universality of 
level fluctuation laws of the spectra of different quantum systems (nuclei, atoms and 
molecules) EI could be understood by the "vacuum" role as a dissipative system [2T| . 
It had already been shown by Maier and Dreizler [22], in a Dirac particle scattering in 
(1+1) dimensions by an electrostatic square well potential, that complex momenta poles 
of the S matrix exhibit the particle-antiparticle content of the Dirac theory and in the 
limit of weak potential strength the poles distribution are similar to the nonrelativistic 
case. 

In summary, we have shown that the zeros of Riemann's zeta function are related 
to the zeros of the Jost function for cutoff potentials in the complex momenta plane in 
a one-to-one correspondence. The energy/width ratios of the large "virtual resonances" 
are associated to the nontrivial zeta zeros and the corresponding widths related to the 
prime sequence. In analogy to the mean field it is expected, by means of a statistical 
hypothesis, that the above relationship would be improved and the distribution of the 
virtual resonances would reflect the chaotic nature of the vacuum. 
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